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Reducing noise in sequence of digitised X-ray images A 
sequence of medical x-ray images is captured over time and 
digitised. The images contain noise which is reduced by filtration. A 
timed noise signal is produced by taking the intensities of a 
particular pixel as they vary with time. A peak (D1) is present and a 
smoothed signal (Ptp) is evaluated using weighting factors 
(W1 ,W2). The weighting factors are applied to the smoothed signal 
at the most recent time. The current noisy signal is modified by the 
probability of the presence or absence of a current noise peak. A 
procedure for evaluating the probability is included. 

Show legal status actions Get Now: Family Legal Status Report 

AT BE OH DE DK ES Fl FR GB GR IE IT LI LU MO NL PT SE 

Show 4 known family members 
Expand full description 

± 1/ DIspositif a rayons X 
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Show all claims 

1. Proc6d6 de fiitrage temporel du bruit dans une image d'une 
sequence d'images num6ris§es en forme de matrices 
bidimensionnelles de pixels, form6es depuis un premier instant (to) 
jusqu'i un dernier instant (t), ces pixels ayant des valours 
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d'intenisit^ correspondant k cheque instant appel6es 6chantillons, 
ce proc6d6 comprenant une evaluation, pour une localisation 
donnee de pixel, d'un 6chantillon du dernier instant temporellement 

filtrd (P^ P), en effectuant une combinaison appel6e premiere 
combinaison de l'6chantillon de Tavaht-dernier instant 

tempbrellement filtr6 (Pc ) affects d'un premier poids (W^) et de 
r6chantillon bruit6 du dernier instant (1^ ^) modul6 par un second 
poids (W ) qui est une fonction d'une valeur de probabilit6 (Pf) 
d'absence de pic de bruit sur cet 6chantillon. 
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(54) Procede de filtrage temporel du bruit dans une image d'une sequence d'images numerisees 
et dispositif mettant en oeuvre ce procede 



(57) L'invention concerne un procede de filtrage 
temporel du bruit dans une image d'une sequence 
d'images numerisees. formdes, d'un premier instant (to) 
^ un dernier Instant (t), de pixels ayant des valeurs d'in- 
tensite correspondant ^ cheque instant appelees 
6chantillons, ce procede comprenant une evaluation, 
pour pixel donne, d'un echantillon du dernier instant 
temporellement filtre (P,p), en effectuant une comblnai- 



son de I'^chantillon filtrd de I'avant-dernier instant (P,« ) 
affects d'un premier poids (W,) et de I'echantillon bruite 
du dernier instant (I,p) module par un second poids (Wg) 
qui est une fonction d'une valeur de probabitite (|5P) 
d'absence de pic de bruit sur cet echantillon. 

Application : Imagerie medicate en fluoroscopie 
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Description 

L'invention concerne un precede de filtrage tempo- 
re! du bruit dans une image d'une sequence d'images 
numerisees, ainsi qu'un dispositif pour mettre en oeuvre s 
ce precede. 

L'invention trouve particulierement son application 
dans le traitement d'images medlcales formees dans un 
mode fluoroscopique de rayons X, au moyen d'un sys- 
teme ou le faisceau de rayons X a une faible intensite, io 
fournissant sous forme de sequence des images nume- 
risees bruitees et peu contrastees qui doivent etre fil- 
trees afin d'en extraire le bruit, sans en'deteriorer les 
details. 

Le filtrage temporel consiste k lisser un signal mo- 
nodlmensionnei appele signal tempore!, constitue par 
les valeurs d'intensite numerisee d'un pixel ayant une 
localisation donnee dans les images de la sequence, en 
fonction du temps. 

Un precede de filtrage temporel est d^ja connu par 20 
la publication de R.E.KALMAN, intitulee "A new ap- 
proach to linear filtering and prediction problems" dans 
"Transactions of the ASME, Journal of Basic Enginee- 
ring, Series 82D, pp. 35-45. 1960. 

Le filtrage de Kalman est defini par une equation 25 
recursive dpnnant I'intensite filtrees d'un pixel courant 
d'une image de la sequence en fonction d'hypetheses 
faites A PRIORI, en fonction de I'intensite du pixel de 
meme localisation dans I'lmage precedente de la se- 
quence et en fonction d'un facteur appele gain de Kal- 30 
man. 

Cette equation peut conduire k plusieurs aigorith- 
mes recursifs. Dans un des algorithmes, le gain de Kal- 
man est une fonction de la difference entre I'intensite 
bruitee observee k un instant donne et I'intensite filtree 55 
a I'instant precedent, pour un pixel de localisation don- 
nee. It en resulte que, s'ii apparait une discontlnuite d'in- 
tensite due k un mouvement ou k un pic de bruit, le si- 
gnal temporel est llsse avant la discontlnuite, mais que, 
par centre, 11 n'est plus filtre apr^s la discontlnuite. Get ^^o 
algerithme n'est done pas approprie au traitement de ce 
genre de discontlnuite d'intensite. 

Un probl^me est que, du k la tr6s faible Intensite du 
faisceau fluoroscopique de rayons X, les Images de la 
sequence sont extremement bruitees, et contiennent ^5 
frequemment des pics de bruit. 

Un autre probldme est que, du au fait que chaque 
image de la sequence est separee de la suivante par 
un petit laps de temps, un evenement important comme 
le deplacement d'un petit outil tel qu'un catheter, peut so 
sun/enir d'une image k Tautre. Le filtrage de I'lmage ou 
apparait le mouvement de ce petit objet ne doit pas de- 
former ou effacer cet objet. 

La presente invention vise k fournir un precede de 
filtrage temporel pour la reduction du bruit, dans les ima- ss 
ges successives d'une sequence d'images 
numerisees : 



qui Iravaille strictement en temps reel, c'est-a-dire 
qui fournisse I'lmage courante filtree des I'appari- 
llon de ladite image courante observee. immedia- 
tement, sans retard d'image, quelle que soit la ca- 
dence des images dans la sequence, 
qui soit capable d'eliminer les pics de bruit : 
qui n'efface pas el ne deferme pas les petits objets 
en mouvement. 

Le prebleme d'elimination des pics de bruit en 
temps reel est resolu au moyen d'un precede de filtrage 
temporel du bruit selon la revendicatlen 1. 

Le prebleme de relimination du bruit sans deformer 
les petits objets en mouvement est resolu par un tel pro- 
cede selon la revendicatlen 5. 

Un dispositif pour mettre en oeuvre ce precede de 
filtrage est recite dans la revendicatlen 7. Ce dispositif 
offre I'avantage d'etre simple a mettre en oeuvre et de 
fournir strictement et exactement en temps reel un fil- 
trage temporel du bruit efficace qui supprime les pics de 
bruit. 

L'invention est decrite ci-apres en detail, en refe- 
rence avec les figures schematiques annexees dent : 

la FIG.1 qui represente un dispositif de prises de 
vues a rayons X : 

la FIG. 2 qui represente une sequence d'images 
numeriques: 

la FIG.3A qui represente un signal tempore! bruite 
cerrespondant a un pixel de localisation donnee 
avec un pic de bruit ; 

ia FIG.3B qui represente un autre signal temporel 
bruite, avec une discontlnuite. due a un 
mouvement : 

la FIG.4A qui represente un signal temporel filtre 
cerrespondant au signal temporel bruite de la FIG. 
3A: 

la FIG. 48 qui represente un signal temporel bruite 
tel que celui de la FIG.4A, d'abord filtre par un filtre 
spatial ; 

les FIGs.SA et 58 qui illustrent un premier et un se- 
cond exemple de filtre spatiaux : 
la FIG. 6 qui montre un exemple possible de fonction 
de modulation de t'echantlllon du present dans la 
fonction de filtrage temporel : 
les FIG. 7 A et 78 qui montrent, sous forme de blocs 
fenctiennels, un dispositif pour mettre en oeuvre le 
precede de filtrage temporel. 

1/ Dispositif a rayons X 

En reference k la FIG.1, dans un exemple d'appll- 
cation du present proc6d6, un systdme de radlegraphie 
numerisee cemprend une source 1 de rayons X : une 
table 2 pour recevoir un patient ; un dispositif intensifi- 
cateur d'image 3 couple k un tube vld6o 4. qui fournit 
des donnees k un systeme 5 de traitement numerique 
d'image comprenant un microprocesseur et approprie k 
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mettre en oeuvre le present precede. Ce microproces- 
seurcomprend plusieurs sorties dont. par exemple, une 
sortie 6 couplee ^ un moniteur 7 pour la visualisation de 
la sequence d'images radiographiques ou des images 
d'intensite. 

La presente invention ne tient pas compte de la me- 
thode par laquelle la sequence d'images numeriques a 
ete obtenue. ni de la nature des objets qu'elle represen- 
te, mais concerne seulement le filtrage de ces images 
danscette sequence pour eliminer le bruit "et les pics de 
bruit. 

11/ Precede de filtrage temporel du bruit 

La presente invention propose un precede de filtra- 
ge temporel du bruit successivement de chaque image 
d'une sequence d'images numeriques bruitees. Ce pro- 
cede opere le filtrage de la derniere image bruitee ob- 
sen/ee. Ce procede est done exactement mene k bien 
en temps reel. 

En reference a la figure 2, ce procede comporte I'ac- 
quisition et la numerisation d'une sequence d'images 
bruitees composee d'images d'intensite du passe Jj^ 
appelees "images causales", formees depuis un pre- 
mier instant jusqu'a un instant j = t-1 , ou t est un instant 
appele present : et d'une image J,^ en cours de filtrage, 
appelee image du present, arrivant k I'instant present j 
= t. Chaque image numerisee bruitee Jj est une matrice 
bidimensionnelle de pixels Aj{x.y) reperes par leurs 
coordonnees x, y el par un niveau d'intensite lj(x, y). 

Comme represente sur les FIG.3A et 3B, un signal 
temporel bruile (It) est forme des differents niveaux d'in- 
tensite l| (X, y) ^ une localisation donnee x.y de pixel Aj 
(x,y), en fonctlon du temps t (lire tau). Ce procede per- 
met d'obtenir un signal temporel filtre P(t), comme re- 
presente sur la FIG.4A. On appelfe "echantiltons" les in- 
tensites qui sont les points des signaux temporels l(r) 
et P(r). Les echantillons du passe ou causals sont notes 
l|^ et I'echantillon du present est note 

Ce procede effectue un filtrage de Techantiilon brui- 
te IjP arrive h I'instant present t, et construit un echan- 
tillon filtrd P,p. Le filtrage temporel est effectue favora- 
blement pour chaque pixel A|(x,y) de la matrice bidimen- 
sionnelle, aux diff^rentes coordonnees x, y, individuel- 
lement. 

La FIG.3A montre un signal temporel bruite I(t) 
ayant un pic de bruit D1 apparaissant d6s I'instant t-1, 
centre sur I'instant t et disparalssant ulterieurement. 

La figure 3B montre un signal temporel I(t) bruite 
ayant une discontinuite d'intensite D2 arrivant k Tinstant 
present t sur I'echantillon I, p. due k un petit mouvement 
local, qui se Iraduit par un front de montee du signal 
temporel I(t) entre I'instant t-1 et I'instant t, entre un si- 
gnal substantiellement bas avant I'instant t-1, et un si- 
gnal qui reste substantiellement haut aprds I'instant t. 

Ce signal temporel montre en outre des petites 
dents de scie dues au bruit, le niveau d'intensite k la 
localisation (x.y) variant constamment. et provoquant un 



aspect de chute de neige dans les images de la sequen- 
ce. Le filtrage permet de supprimer cet aspect en lissant 
a une valeur moyenne le signal temporel I(t) relatif a un 
pixel donne. Ce filtrage a en outre pour but -d'eliminer 

5 les pics de bruit tels que 01 . 

Un pic de bruit tel que 01 doit etre distingue d'une 
discontinuite d'intensite telle que 02 pour pouvoir etre 
elimine. A cet effet. le present procede prend en compte 
que le pic de bruit 01 ne correspond pas a un pheno- 

10 mene purement spatial comme cela est le cas de la dis- 
continuite 02, mais correspond en outre a un phenome- 
ne temporel. 

L'echantillon filtre P,^ pour reconstruire I'echantillon 
bruile du present 1,^ s'exprime par la relation (1) ; 

IS 



20 Selon cette relation ( 1 ) I'echantillon filtre P^^, ou si- 
gnal temporel filtre k I'inslanl present t, est une combi- 
naison de I'echantillon causal filtre P^^ evalue a I'avant- 
dernier instant, appele dernier echanlillon causal filtre. 
et de I'echantillon bruite 1^^ observe a I'instant present 
25 I Dans la formule (1 ) du signal filtre PjP. on affecte aux 
echantillons P^^^ et 1,^ respectivement un premier poids 
et un second poids Wg. 

Le denomlnateur de la formule (1 ) est un facteur de 
normalisation dont la presence est fondee sur le fait que 
30 la somme des poids que Ton applique aux echantillons 
doit etre egale a 1 pour que la moyenne du signal filtre 
PjP soil egale k la moyenne du signal bruite I, p. 

Dans la formule (1) I'echantillon bruite du present 
l^^ est module par une fonction d'une valeur de proba- 
35 bilite (5^ d'absence de pic de bruit a I'instant present I 
A cet effet, une detection de pic de bruit est appliquee 
au signal temporel a I'instant present t, pour evaluer la 
valeur de probabilile (3^ d'absence de pic de bruit, et le 
poids affecte k I'echantillon bruite du present 
^0 dans la combinaison (1 ) est de preference une fonction 
de cette valeur de probability p^. Dans un exemple ap- 
proprie, ce second poids est tel que ; 

4S W2 = |i = f (2) 

En reference k la FIG.7A, qui illustre le present pro- 
cede au moyen de blocs fonctionnels, la detection de 
pic de bruit sur I'echantillon bruite du present l,P 
50 comprend : 

un passage 1 01 d'un filtre spatial G sur I'echantillon 
bruite du present Ij^. pour evaluer un echanlillon lis- 
se du present Gj^, 
55 una evaluation 1 03 d'une premiere difference A1 ef- 
fecluee entre rechantillon bruite du present I,p el 
rechanlillon lisse du present G,p : 
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A1 =1,^0,'' (3) 

une evaluation 104 d'une seconde difference A2 ef- 
-fectuee entre I'echantillon bruite du present et le 
dernier echantillon causal bruite ; 

^2 = 1,^ - ij:, (4) 

une connparalson 107 des signes de deux differen- 
ces et S2. 

une evaluation 1 1 0 d'une variable superieure ou 
egale ^ zero telle que Z,p > 0. L'evaluation de la 
variable Z, ^ est une nnesure de presomption de pre- 
sence de pic de bruit a I'instant present t. En effet, 
un pic de bruit est un phenomene transitoire a la 
fois dans le domaine spatial et dans le domaine 
temporel. 

L'evaluation 103 de la prenniere difference A1 entre 
I'echantillon bruite 1,^ ^ I'instant t et I'echantillon filtre 
spatialement G,p relatif au meme instant permet de met- 
tre en evidence I'aspect transitoire dans le donnaine spa- 
tial. 

Pour realiser le filtrage spatial de I'echantillon bruite 
a I'instant du present t, toute fonction G de filtrage spatial 
qui supprime les pics locaux de signal, ou les extrema 
locaux de signal, est appropriee a etre utiiisee. Ce filtra- 
ge spatial peut etre du type passe-bas, ou bien de pre- 
ference du type median hybride tel que decrit par exem- 
ple dans la publication intitulee "A new Class of Detail 
Preserving Filters for Image Processing" par Ari NIEIVII- 
NEN et alii, dans IEEE Transactions on Patterns Analy- 
sis and Machine Intelligence, vol. PAMI.9. N**1, Jan. 
1987. De preference, selon ['invention, un des filtres 
spatiaux 101 tels que decrits ci-apr6s, sera applique au 
signal bruite l^^. 

En reference k la FIG.5A, dans un premier exemple 
de filtre spatial 1 01 , on fait passer sur I'image Jj a traiter 
et au pixel A^^x, y), un filtre median de quelques pixels, 
par exemple 3 pixels. Dans une petite fenetre U, (x, y). 
ce filtre 101 centrd sur le pixel courant A^ (x. y) calcule 
I'intenslte des 3 pixels adjacents, sur une direction don- 
nee, par exemple horizontale, ordonne leurs valeurs 
d'Intensite de la plus petite k la plus grande intensite et 
retient comme intensite filtree, Tintensite autre que la 
plus petite el la plus grande. appel^e intensity mediane. 
Cette intensite m^dlane GjP est affect6e au pixel cou- 
rant At (x. y). 

En reference h la FIG.5B dans un second exemple 
de filtre spatial 1 01 . on fait passer sur I'image ^ traiter 
el au pixel A, (x, y), un filtre median d6limite par une 
fenetre {x, y) centree en x. y ce filtre ayant plusieurs 
orientations possibles par exemple une orientation ho- 
rizontale appelee H. une orientation verticale appelee 
V. el des orientations selon une premiere diagonale ap- 



pelee XI, et selon une seconde diagonale appelee X2. 

La valeur mediane est d'abord determinee surcha- 
que orientation. Puis deux nouvelles valeurs medianes 
sont determinees en prenant le resultat des deux orien- 
s tations H et V avec I'intensite du point central A, (x, y) 
d'une part, et le resultat des deux orientations XI et X2 
d'autre part avec I'intensite de point central A^ (x, y). En- 
fin, une dernlere valeur mediane est determinee en pre- 
nant ces deux valeurs nouvelles avec I'intensite de point 
10 central A, (x. y). Cette dernlere valeur mediane G^^ est 
affectee au point central A, (x. y). 

L'evaluation 104 de la seconde difference :12 entre 
les deux echantlllons consecutifs dans le temps formes 
par I'echantillon I, ^ bruite a I'instant t et Techantillon brul- 
^5 te 1,^^ \ I'instant precedent t-1 , permet de mettre en evi- 
dence I'aspect transitoire dans le domaine temporel. 

La presence d'un phenomene transitoire ^ la fois 
dans le domaine spatial et dans le domaine temporel. 
c'est-^-dire la presence d'un pic de bruit, est detectee 
20 lorsque les differences ^^ et A2 evaluees sont coheren- 
tes et monlrent une variation d'Intensite dans le meme 
sens, ce qui est concretise par une identite des signes 
de ces deux differences A1 . ^2. 

La presence de pic de bruit est presumee nulle lors- 
25 que les signes des differences :\^ . A2 sont differenls : 

Z,^ - 0 (5) 

^0 La presence de pic de bruit est presumee averee 
lorsque les signes des differences ai, a2 sont les 
memes : 

35 Z^%0 (6) 

La presomption de presence de pic de bruit est 
alors mesuree selon I'lnventlon par la relation : 

40 

Z, =f\/1in(Ul|. |A2|) (7) 

A cet effet. le proced6 comprend revaluation t08; 1 09 
des valeurs absolues des differences M et ii2 

45 

lAlNl't^'-Gtl . (8) 

50 U2i=|i^i;:,| (9) 

et la mesure Zj^ de presomption de presence de pic de 
bruit est prise comme la plus petite valeur enlre les deux 
valeurs absolues des differences A1 . ^2. 

Le precede comprend alors revaluation 110 d'une 
probablllte d'absence de pic. Cette probabilite ji^ est 
evaluee comme une fonction F.^ non strlclement d6- 



BNSOOCtO: <EP oagiaggAI I > 



4 



EP 0 821 322 A1 



' croissante de la variable Z,'' selon la relation : 



:0.1 



(13) 
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(10) 



Pour eviter que le precede de filtrage temporel, in- 
cluant ces etapes de detection de pic de bruit et d'eva- 
luation de probabilite ne donne des resultats aberrant 
dans les cas ou en meme temps : 

il existe une discontinuite d'intensite reelle dans le 
passe, ce qui implique que le precede ne peut pas 
s'appuyer sur les echantillons causals precedant la 
discontinuite. 

il existe un pic de bruit a Tinstant present, ce qui 
innplique que le precede ne.peut pas s'appuyer sur . 
I'echantillon present, alors on interdit ^ la fonction 
decroissante de prendre la valeur zero. Cela si- 
gnifie que, selon le precede, la probabilite d'ab- 
sence de pic a I'instant courant t n'est jamais nulle 



(11) 



En reference a la FIG. 6, une courbe representative 
possible de la probabilite en fonction de la va- 
riable montre a titre d'exemple que : 

la probabilite d'absence de pic vaut 1 lers- 
que la valeur Z^^ de presomption de presence 
de pic vaut 0, 

la probabilite [5^ d*absence de pic reste proche 
de 1 pour les valeurs Z^^ de presomption de 
presence de pic inferieures ^ 1 , 
la probabilite p'^ d'absence de pic tend vers une 
valeur pour les valeurs Z^p de presomp- 
tion de presence de pic superieures a 2. 

Dans la formule (1), le poids de I'echantillon du 
present I^p est pris egal k la probabilite ji*^ d'absence de 
pic. Lorsque la probabilite d'absence de pic est 1, il y 
certitude d'absence de pic, et la confiance qui doit etre 
accordee ^ I'echantillon du present I^p est grande. Pour 
traduire ce fait, le poids W2 associe k I'echantillon du 
present IjP a la valeur 1 cbmme la probabilite soit : 



W2 = 13,^ = 1 



(12) 



Lorsque la probabilite d'absence de pic (i^ est de I'ordre 
de I'ordre de li^iN- c'est-^-dire a une valeur faible, il y a 
quasi-certitude de presence de pic et la confiance ac- 
cordee a I'echantillen du present 1,^ est faible. Pour tra- 
duire ce fait, le poids associe ^ I'echantillon du present 
prend une valeur faible, par exemple : 
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Entre ces deux valeurs. I'echantillon du present Ij^* est 
module, dans la formule (1). par la valeur de la proba- 
bilite p*^ calculee comme decrit plus haut en fonction des 
differences :^^ et A2, 

II est souhaitable que la courbe (3'j' en fonction de 
la variable Z, ^ ne depende pas du bruit. Pour obtenir ce 
resultat. chacune des differences A1 et :\2 est nermali- 
see par I'ecart type du bruit Og, qui est une notion con- 
nue de I'homme du metier. 

On rappelle que I'ecart type du bruit est la racine 
carree de la variance du bruit a| relative ^ un echantillon 
donne. La variance du bruit est la meyenne des ecarts 
au carre du bruit dans la region consideree de Timage 
Jj par rapport a la meyenne locale d'intensite dans cette 
region. 

Pour chaque difference . A2, la variance du bruit 
est approximativement 2a| . II en resulte que I'ecart type 
du bruit est approximativement Og V2. 

Les differences utilisees sent dene de preference 
les differences nermalisees A1 et A2 dennees par les 
relations : 



^2 : 



0bV2 



(14) 



f15^ 



D'eu il resulte que : 



: Min ( 



IA1I lA2t 



■'8 



V2' 



a^J~2 



(16) 



Le precede decrit plus haut permet de distinguer un 
pic de bruit comme montre sur la FIG.3A d'une discon- 
tinuite comme montree sur la FIG. SB. permet de pren- 
dre en compte un tel pic de bruit, et d'eliminer ce pic de 
bruit comme montre sur la FIG.4A. 

Dans ce precede, le dernier echantillon causal filtre 
doit dtre calcule par un filtrage temporel qui peut ^tre un 
precede recursif ou non recursif connu de Thomme du 
metier De preference, le dernier 6chanti!lon causal est 
calcule par une seconde combinaison telle que : 



^ c-i 



j=t-i 
j -to 



(17) 
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ou Ij^ represente les echantillons causals bruites, pre- 
leves dans les images de la sequence, ^ une localisa- 
tion A,(x,y) donnee de pixel, du premier instant to, jus- 
qu'^ Tavanl-demier instant t-1 , et ou bj^ represente des 
poids associes respectivement ^ chaque echanlillon de 
m6me indlce. La seconde comblnaison est normalisee 
par la somme des poids bj^. 

Dans un exemple de mise en oeuvre du procede, 
les poids bj^ peuvent etre des constantes affectees res- 
pectivement aux echantillons causals, determinees em- 
piriquement par I'homme du metier. Dans ce cas, ce- 
pendant, le probl6me du mouvement des tres petits ob- 
jets n'est pas completement resolu. 

De preference, un poids bf^^ d'un echantillon donne 
causal Ik^ est le produit des probabilites d'absence de 
discontinulte a9, effectue entre cet echantillon causal 
donne 1,^^ et I'echantillon bruite du present 1,^. 

Par discontlnuite, on entend une discontinulte d'in- 
tensite entre deux echantillons consecutifs de la suite 
d'echantillon causals, due ^ un mouvement. et done liee 
a un phenomene spatial, comme par exemple montre 
sur la FIG.3B .en 2D. Ainsi : 



j-t-i 

W3 = E 

j-co 



(21) 



D'une manl^re preferee, on pose, k litre de simpli- 
fication, I'expression : 



10 



IS 



20 



(22) 



L'expresslon de la formulation (1) de la premiere 
combinaison pour obtenir le filtrage des pics de bruit uti- 
lise alors de maniere avantageuse I'expression d'un troi- 
sieme poids W,p donnee par : 



(23) 



3-C 

= n 



(18) 



25 de telle maniere que revaluation de I'echantillon filtre du 
present soil donnee par la relation : 



A titre d'exemple, la formulation du poids causal lie 



S I'echantillon causal bruite l,p, est : 



b,., ^ ot, 



(19) 



30 



qui est eauivalente k la formulation : 



qui est la probabilite d'absence de discontinulte entre 55 
les echantillons et 1,^. Et la formulation du poids cau- pu bien 



sal b,^2 ''^^ I'echantillon causal bruite \^2^^X 



.COG 

b,.2 = a, X a,.^ 



(20) 



40 



(24a). 



(24b) 



(1) 



etc. 

II resulte de cette formulation, que si une disconti- 
nuity d'intensite telle que D2 inlervient entre I'instant ou 
arrive I'echantillon 1,^^ et I'instant ou arrive Techantillon 
du present I, p. alors ledit echantillon est peu pris en 
compte dans la seconde combinaison car son poids est 
faible. Seuls les Echantillons pour lesquels il y a conti- 
nuity d'intensite jusqu'^ Techantillon du present sont for- 
tement pris en compte dans la seconde combinaison de 
la relation (17). 

Selon le procEdE de filtrage temporel du bruit et des 
pics de bruit, dans la formule (1 ). le dernier echantillon 
causal filtre est dvaluE selon la formule (1 7) et le poids 
qui lui est assocle est donnes par ladite somme des 
poids de la seconde combinaison : . 



ou encore equivalente k I'expression : 



(25) 



qui montre que dans la premiere combinaison (1) 
r^chantiilon du present 1,^ est moduli par un coefficient 
de probabilite (i^= W2 d'absence de pic de bruit et que 
I'avant-dernier Echantillon causal P^^ a pour poids 
la somme des poids 1/K,f^ de la seconde combinaison 
(17) qui permet la determination dudit Echantillon. 

En utilisanl la mEthode preferentielle qui vient d'Etre 
dEcrite pour la dEtermination de Tavant -dernier echan- 
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'lillon causal filtre et de son poids W1 Je probleme de 
relimination du bruit sans deformer ou effacer les petits 
objets en mouvennent provoquant des disconttnuites 
spatiales telles que D2 de la F1G.3B. est completement 
resolu. 5 

Les premiere, seconde et troisieme combinaisons 
peuvent etre respectivement linealres. 

Ill/ Dispositif pour mettre en oeuvre le precede de filtra- 

La formule (24b) est aisement evaluable par des 
moyens de calcul simple, comme il est decrit ci-apres. 

En reference aux FIG.7A et 7B, un dispositif simple 
pour mettre en oeuvre le precede de filtrage propose 75 
est decrit sous fornae de blocs fonctlonnels. 

Ce dispositif comprend d'abord des moyens non re- 
presentes, pour extraire les echantillons bruites relatifs 
^ un pixel Aj(x.y) de localisation donnee x,y dans les 
images numeriques, du premier instant to jusqu'a Tins- 20 
tant du present t. 

En reference a la FIG.7B, ce dispositif comprend 
ensuite un premier ensemble 1 90, que Ton peut appeler 
ensemble de filtrage temporal causal, pour effectuer 
revaluation, a partir des echantillons causais bruites 2S 
preleves dans le temps de to a t-1. a la localisation don- 
nee x.y, du dernier echantillon causal filtre P^^ par un 
filtrage lineaire, avec des poids associes aux echan- 
tillons qui sont, soit constants, soit de preference fonc- 
tions des coefficients a9 de probabilite d'absence de 30 
discontinuite comme decrit plus hauL 

Des fillres recursifs ou non recursifs. appropries ^ 
realiser un filtrage lineaire. sont connus de Thomme du 
metier et ne sont pas decrits plus avant. Un tel filtre li- 
neaire 1 90 pour effectuer ce filtrage causal a une entree 35 
189 pour les echantillons causats bruites l|C et une en- 
tree 188 pour les poids bj^ constants ou sous la forme 
preferee de produits des a9et fournit ^ sa sortie 204 le 
dernier echantillon causal riltre P^^ , et a sa sortie 205, 
dans la forme preferee, la somme des poids 40 

sous la forme inverse qui est posee egale ^ K^*^^ . 

Ce dispositif comprend ensuite des moyens ctassi- 
ques non representes, car ^ la portee de tout homme 
du metier, pour evaluer Tecart type du bruit Oq v2 ^ partir 50 
de la variance du bruit a| des echantillons des differen- 
ces ^^ ,A2. 

En reference ^ la FIG.7A. le dispositif de filtrage 
temporel pour diminuer les pics de bruit comprend es- 
sentiellement un ensemble 200 incluant des entrees qui ss 
sont : 

une entree 99 pour le dernier echantillon bruite cau- 



une entree 100 pour I'echantillon bruite du present 

1,^ 

une entree 204 pour le dernier echantillon causal 
filtre Pj5^ provenant du filtre lineaire causal 190, 
une entree 205 pour la vafeur K^^ provenant du fil- 
tre lineaire causal 190, 

une entree 206 pour les valeurs evaluees de I'ecart 
type du bruit ag, 

une sortie 301 pour I'echantillon filtre du present P,p 
ou le bruit et les pics de bruit sont elimines de fagon 
extremement rigoureuse. 

Dans I'ensemble 200, se trouvent : 

un premier sous-ensemble 21 0 pour evaluer la pro- 
babilite [J*^ d'absence de pics de bruiL 
et un deuxieme sous-ensemble 220 pour evaluer 
finalement Techantillon temporellement filtre du 
present Pj^. 

Le premier sous-ensemble re^oit I'echantillon bruite 
du present I[P et le filtre spatialement par le filtre spatial 
G, en 101, pour fournir I'echantillon lissedu present GjP 
Ce premier sous-ensemble regioit aussi le dernier 
echantillon causal filtre t(?^ et la valeur evaluee de I'ecart 
type du bruit ag. 

L'echantillon bruite du present l^^ et I'echantillon lis- 
se du present sont portes dans un premier addition-" 
neur 103 qui fournit la premiere difference A1 ; I'echan- 
tillon bruite du present \f et le dernier echantillon causal 
bruite sont portes dans un deuxieme additionneur 
104 qui fournit !a deuxieme difference Les differen- 
ces A1 et d2 sont entrees dans des Tables (LU.T) res- 
pectivement 105 et 106, de meme que la valeur de 
I'ecart du bruit ag pour fournir lesdites differences nor- 
malisees par I'ecart du bruit selon les relations respec- 
tives (14) et (15). 

Le signe des differences normalisees est compare 
dans le comparateur de signe 107 qui fournit un signal 
0 ou 1 pour imposer que la valeur de la mesure de pre- 
somption de presence de pic de bruit soit evaluee : 

0 si les signes sont diff^rents, 

et 9i 0 si les signes sont identiques. 

Les valeurs des differences A1. A2 normalisees 
sont portees respectivement dans des blocs 108, 109 
d'^valuation des valeurs absolues : ^ la suite de quoi. 
la probability est calculee dans la Table (LUT) 110, 
selon la commande du comparateur 107, comme urie 
fonction F^^ non strictement decroissante de la valeur 
de presomption 2,p selon la relation (10). 

Le second sous-ensemble 220 re^oit I'echantillon 
bruite du present l^^, le dernier Echantillon causal filtr6 
P,p^ . et rinverse de la somme des poids causais not6e 
K(^^ comme ddj^ decrit. Ce second sous-ensemble 220 
comprend une Table 111 (L.U.T) ou sont portees la va- 
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leur de la probabilite d'absence de pic 15"^, et la valeur 
de I'inverse de la somme des poids causals K^^ , pour 
fournir le poids W,p evalue selon la relation (23). 

Ce second sous-ensemble comprend ensuite des 
additionneurs 112, 114 et un multiplieur 113 pour eva- s 
luer I'echantillon filtre du present P^^ en ulilisant la re- 
lation (24b). A cet eWel I'additionneur 1 1 2 re(?oit I'echan- 
tillon bruite du present \^^ ei le dernier echantillon causal 
filtre P^^ et fournit la difference - P,p^ qui est multi- 
pliee par le poids W^p dans le multiplieur 113 dont la 
sortie est ajoutee, par le moyen de I'additionneur 1 1 4 au 
dernier echantillon causal filtre Pj^^ . II en resulte la va- 
leur de I'echantillon filtre du present P,'' disponible a la 
sortie 301 du sous-ensennble 220, confornnement a la 
formule (24b). is 



Revendications 

1 . Precede de filtrage tennporel du bruit dans une ima- 20 
ge d'une sequence d' images numerisees en forme 

de matrices bidimensionnelles de pixels, formees 
depuis un premier instant (to) jusqu'a un dernier ins- 
tant (t), ces pixels ayant des valeurs d'intensite cor- 
respondant k chaque instant appelees echantlllons, 2$ 
ce procede comprenant une evaluation, poui* une 
localisation donnee de pixel d'un echantillon du 
dernier instant temporellement filtre (P^P). en effec- 
tuant une combinaison appelee premiere combinai- 
son de i'echantillon de I'avant -dernier instant tern- 30 
porellement filtre (P,^, ) affecte d*un premier poids 
(Wi) et de I'echantillon bruite du dernier instant (I,P) 
module par un second poids (Wg) qui est une fonc- 
tion d'une valeur de probabilite (|5^) d'absence de 
pic de bruit sur cet echantillon. 35 

2. Procede selon la revendication 1 , comprenant une 
evaluation d'une mesure de presomption de pre- 
sence de pic de bruit (Z,^), superieure ou egale a 

0, comme la valeur m in imale entre une premie re va- -^o 
leur absolue de la difference entre I'echantillon brui- 
te du dernier instant (I,P) et cet echantillon spatia- 
lement lisse (G,P), et une deuxieme valeur absolue 
de la difference entre les echantillons bruites des 
dernier (1^^) et avant-dernier instants (1^^^ ). et une ^5 
determination du second poids (Wg) par evaluation 
de ladlte fonction de la valeur de probabilite (|i^) 
d'absence de pic de bruit sur I'echantillon du dernier 
instant comme une fonction (F^^) non strictement 
decroissante, superieure ^ 0 et inf^rieure ou egale so 
k 1, de la mesure de presomption de presence de 
pic de bruit (Z^P). 

3. Procede selon la revendication 2, selon lequel cha- 
cune des premiere et seconde valeurs absolues 55 
des differences est normalisee par un ecart type du 
bruit evalue sur les Echantillons pris en compte 
dans lesdites differences. 



4. Procede selon I'une des revendications 1 a 3. com- 
prenant une evaluation de I'echantillon temporelle- 
ment filtre de I'avant-dernier instant (P^^^ ) en effec- 
tuant une combinaison appelee seconde combinai- 
son des echantillons bruites (Ij^) du premier instant 
a I'avant-dernier instant de la sequence, affectes 
respectivement de poids appeles poids causaux 
(bj^). et ladite seconde combinaison etant normali- 
see par la somme desdits poids causaux. 

5. Procede selon la revendication 4, comprenant une 
evaluation du poids (W,) pour appliquer a I'echan- 
tillon (P^^ ) temporellement filtre de I'avant-dernier 
instant dans ladite premiere combinaison comme la 
valeur de la somme des poids causaux 

j-t-i 

( E = wj 

j -co 

de la seconde combinaison. et une normalisation 
de ladite premiere combinaison par la somme des 
premier et second poids (W-i+Wg). 

6. Procede selon la revendication 5, comprenant une 
evaluation d'un troisieme poids (W,^) donnee par 
un quotient du second poids (Wg) par la somme des 
premier et second poids (W^+Wg), et une evalua-' 
tion de Techantillon temporellement filtre du dernier 
instant (P,^) en effectuant une combinaison appe- 
lee troisieme combinaison de I'echantillon (P,^^) 
temporellement filtre de I'avant-dernier instant et 
d'une difference de I'echantillon bruite (l,^) du der- 
nier instant et de I'echantillon {P^^ ) temporellement 
filtre de I'avant-dernier instant, difference mod u lee 
par ledit troisieme poids (Wt^). 

7. Dispositif pour mettre en oeuvre un procede selon 
I'une des revendications 1 k 6, comprenant un sys- 
teme pour fournir une suite d'echantillons bruites (Ij) 
relatifs a une localisation donnee de pixels dans des 
images successives d'une sequence d'images 
ayant la forme de matrices bidimensionnelles de 
pixels, depuis un premier instant (to) jusqu*^ un der- 
nier instant (t), un premier ensemble ( 1 90) de calcul 
pour fournir un echantillon filtre temporellement 
(Pt?i) cie I'avant-dernier instant (t-1), et un second 
ensemble (200) de calcul comprenant un premier 
sous-ensemble (210) incluant des modules de cal- 
cul pour evaluer une fonction d'une valeur de pro- 
bability (p*^) d'absence de pic de bruit sur I'echan- 
tillon bruite du dernier instant (1,^), et un second 
sous-ensemble (220) incluant des modules de cal- 
cul pour ^valuer une premiere combinaison de 
I'echantillon temporellement filtre (P,f^ ) de I'avant- 
dernier instant affecte d'un premier poids (W, ) el de 
I'echantillon bruite (1,^) du dernier instant affecte 
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d'un deuxi^me poids (W2) determine par la fonctton 
de la valeur de probabiltte ), cette premiere com- 
binaison etant normalisee par la somme du premier 
et du deuxieme poids. 

5 

8. Dispositif selon la revendication 7, comprenant 
dans le premier ensemble de calcul (190), des 
moyens pour evaluer rechantillon temporellement 
filtre (P,p^ ) de Tavant-dernier instant comme une se- 
conde combinaison des echantillons' bruites (IjC) w 
anterieurs au dernier instant affectes de poids res- 
pectifs (bjC) appeles poids causaux, et normalisee 
par la somme desdits poids causaux, et compre- 
nant dans te deuxieme ensemble de calcul (200), 
des moyens pour evaluer le premier poids (W^ ) af- 
fecte ^ I'echantillon temporellement filtre {P^^ ) de 
I'avant-demier instant comme la somme desdits 
poids causaux 



20 
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Dispositif pour le traitement d'images medicales 25 
comprenant un systeme pour fournir des donnees 
d'une image sous forme de matrice bidimensionnel- 
le de pixels ayant des valeurs d'intensite numeri- 
sees, un systeme d'affichage des donnees, un sys- 
teme de traitement d'image ayant acces aux don- 30 
nees d'image et au systeme d'affichage, et un dis- 
positif selon Tune des revendications 7 ou 8. 
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